Mass Outflow in Active Galactic Nuclei: New Perspectives 

ASP Conference Series, Vol. TBD, 2001 

D.M. Crenshaw, S.B. Kraemer, and I.M. George 



A physically consistent model for X-ray emission by 
Seyfert 2 galaxies demonstrated on NGC 1068 

E. Behar, A. Kinkhabwala, M. Sako, F. Paerels, S.M. Kahn 
Columbia University, 550 W 120th St., New York, NY 10027 

A.C. Brinkman, J. Kaastra, R. van der Meer 

SRON, Sorbonnelaan 2, 3548 CA, Utrecht, The Netherlands 

Abstract. Preliminary analysis of the X-ray spectrum of NGC 1068 
obtained by the RGS spectrometer on board XMM-Newton is presented. 
A physically consistent model is developed in order to quantitatively de- 
scribe the reprocessing of the central AGN continuum source into the 
discrete X-ray emission observed in Seyfert 2 galaxies. All the important 
atomic processes are taken into account, including photoexcitation, which 
has been neglected in some previous models. The model fits the high res- 
olution NGC 1068 data very well, which implies that the contribution of 
hot collisional gas to the X-ray spectrum of NGC 1068 is negligible. 



1. Introduction 

Recent high-resolution X-ray observations of Seyfert 2 galaxies with the grating 
spectrometers on board XMM-Newton and Chandra reveal an emission spectrum 
rich in discrete features. This finding supports the unified AGN model, accord- 
ing to which the continuum nuclear source in Seyfert 2 galaxies is completely 
blocked from our line of sight and we can see only the reprocessed light. To date, 
discrete emission has been detected in Markarian 3 (Sako et al. 2000), NGC 4151 
(Ogle et al. 2000), and the Circinus Galaxy (Sambruna et al. 2001). NGC 1068, 
which is the brightest Seyfert 2 and therefore best suited for detailed spectro- 
scopic studies, has been observed by all X-ray grating spectrometers: LETGS 
(Brinkman et al. 2001), HETGS (Ogle et al. 2001), and RGS (Kinkhabwala et 
al. 2001 and the present work). The conspicuous radiative recombination con- 
tinuum (RRC) features observed in all Seyfert 2 galaxies verify the photoionized 
nature and low-temperature {kTf> ~ few eV) of the X-ray emitting plasma. 

Spectra emitted by photoionized plasmas have been occasionally modeled 
assuming that the spectral lines and RRC's are solely produced by recombination 
and subsequent radiative cascades (e.g., Liedahl et al. 1990, Sako et al. 1999, 
Porquet & Dubau 2001). When comparing these pure-recombination models to 
the Seyfert 2 X-ray spectra, significant residuals are found, which are mostly in 
lines with high oscillator strengths that are also among the strongest lines in 
hot collisional plasmas. Sako et al. (2000) were the first to suggest that these 
lines arise from photoexcitation (resonant scattering), while Ogle et al. (2000, 
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2001) interpret them as evidence for hot (^-^keV) plasma confining the much 
colder photoionized gas. Neither explanation has been tested with a detailed 
model yet. It is only natural that in a strong radiation field environment, the 
same flux that ionizes the atoms (bound-free transitions) will also excite them 
(bound-bound). In fact, emission by recombination and subsequent cascades 
without photoexcitation would be possible only if the plasma is in a transient 
phase of cooling and recombining. However, as a general rule in ionization bal- 
ance, both of these processes contribute appreciably to the emitted spectrum. 
The relative importance of photoexcitation, recombination, and coUisional ex- 
citation can be assessed only through careful modeling of the various line (and 
RRC) intensities, while correctly accounting for the atomic state kinetics. The 
model suggested in this paper is similar in spirit to the models that have been 
employed to interpret the spectra emitted by laboratory laser-produced plasma 
(Doron et al. 1998) and stellar coronae (Behar, Cottam, & Kahn 2001). Only 
here, we add the transitions induced by the radiation field, namely photoioniza- 
tion and photoexcitation. We focus on the unique spectroscopic signatures of 
photoexcitation and point out their importance in the NGC 1068 spectrum ob- 
tained with the Reflection Grating Spectrometer (RGS) on board XMM- Newton. 
The full analysis of the NGC 1068 RGS observation will be presented elsewhere 
(Kinkhabwala et al. 2001). 



2. The Atomic-State Kinetic Model 

For each ion, we use a steady-state model, which includes all of the excited 
levels of the emitting ion in addition to the ground level of the next charge 
state. For the simple case of H-like ions, that would be the nl configurations in 
addition to the bare nuclei state. The atomic processes included in the model 
are radiative decays, radiative recombination (RR), photoexcitation (PE), and 
photoionization (PI). Additional processes, such as collisional excitations, de- 
excitations, or dielectronic recombination processes can be easily incorporated, 
but their effect in typical Seyfert 2 plasmas is negligibly small. Self absorption 
of the reprocessed light can be included, but has been found not to be important 
for NGC 1068. The electron energy distribution is assumed to be Maxwellian, 
corresponding to an electron temperature Tg. The set of rate equations for the 
density in cm~^ of an ion with positive charge g in a level i can be written 
as: 

j>i k>i j<i 

-nt' f E + E + E = (1) 

\j<i j>i k>i / 

The electron density is denoted by He and A represents the Einstein coeffi- 
cient for spontaneous radiative emission. a^^(Te) is the rate coefficient for RR. 
The rates R^^ for PE are calculated from the illuminating photon flux F{E), 
which is a function of the photon energy E and is expressed in photons sec~^ 
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The cross section for PE from level j to level i is: 



,2 




nieC 



km 



(3) 



Here, e and rUe are the electron charge and mass, c is the speed of light, fji is the 
oscillator strength of the line, and (t){E) is the line profile. Analogously, the rates 
R^^ for PI are calculated using the PI cross section and the same flux F(E). 
F(E) is taken to be a power law [i^o(-^)] absorbed along the ionization cone. The 
optical depth t{E) is calculated taking both line (PE) and edge (PI) absorption 
into account. Since with the RGS, we measure the spectrum integrated over the 
entire ionization cone of NGC 1068, we use for Eq. (1) the effective flux : 



For spatially resolved spectra obtained with the Chandra gratings, a different 
flux needs to be used for each region. We will present a more detailed description 
of the model in Kinkhabwala et al. (2001). 

Level-by-level calculations are carried out for each ion that appears in the 
spectrum. High lying levels with n-values up to 7 are included explicitly in the 
model. The contribution of higher levels, both via RR and via PE towards these 
levels, is included by extrapolation. The set of Eqs. (1) is solved by normaliz- 
ing the sum of all level populations to unity. The power law normalization (or, 
alternatively, the electron density) is tuned to give the appropriate ionization 
balance between the two adjacent ions, which reproduces the observed RRC to 
line ratios. The line intensities in number of photons emitted per second and 
per cm^ are finally given by n^''Aij. The elementary atomic quantities used in 
this work, namely level energies, radiative decay rates, and photoionization cross 
sections, were obtained by means of the multi-configuration, relativistic HUL- 
LAC (Hebrew University Lawrence Livermore Atomic Code) computer package 
developed by Bar-Shalom et al. (2001). 

3. Signature of Photoexcitation in the High Rydberg Series Lines 

Looking at a Rydberg series of lines, for instance the Lyman series, the ratios of 
high-n lines to that of, say, Lya provide a probe for PE. At high ionic column 
density, beyond that in which the absorption of Lya is saturated, the high-order 
lines can still be enhanced by PE resulting in the increase of their intensity ratio 
to Lya. This effect is illustrated in the left panel of Fig. 1 for C"*"^. At the 
point where Lya saturates (~ lO^'^^cm^^), the ratios rise, eventually reaching 
values of up to a few times higher than the low column-density values. As the 
column density increases further, the other members of the series progressively 
saturate as well. The actual column densities for which all of these saturations 
occur depend on (p{E) (essentially the turbulent velocity, taken here to be 800 



F{E) = Fo{E) 



1 - exp[-T{E)] 



(4) 
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Figure 1. Left- Calculated intensity ratios of the Lyman series lines 
and RRC (circles) of C"*"^ to the Lya intensity. The NGC 1068 values 
for these ratios, as obtained from the RGS spectrum, are also plotted 
(asterisks), as well as the values for collisional gas at 100 cV (hor- 
izontal bars). Right- The NGC 1068 spectrum (dotted line) obtained 
by the RGS spectrometer on board XMM-Newton in the region of the 
C+^ lines and RRC compared with the present model (solid line) at 
10^^ cm-2. 

km/sec), but the intensity ratio trend seen in Fig. 1 is generic to the Seyfert 2 
scenario, in which photons from a power law source are reprocessed by means 
of PE and PI. Next to the theoretical ratio curves in Fig. 1, still in the left 
panel, the ratios measured for NGC 1068 are shown (asterisks). These ratios 
are clearly in the high ionic column density regime at about 10^^ cm~^. On 
the right hand side of the left panel in Fig. 1, the same ratios are indicated 
for collisional conditions where C+'' forms (fcTg = 100 eV). Clearly, these are 
lower than the ratios for photoionized conditions, even at their low column- 
density limit. Hence, an additional collisional component will only worsen the 
agreement of the model with the observed ratios. 

Indeed, the current, rather simple model for 10^^ cm~^ fits the measured 
spectrum of NGC 1068 fairly well, as can be readily seen in the right panel of 
Fig. 1. The small remaining discrepancies could be due to the approximate 
line profiles used here. The complete analysis of the NGC 1068 spectrum, which 
includes all of the observed ions as well as the corrections for the line profiles and 
outflow velocity shifts will be presented in Kinkhabwala et al. (2001). Here, the 
slight velocity shifts actually come in handy and assist the comparison between 
the data and the model. 



4. He-like species 

The He-like species feature, in addition to the leading Is - np Rydberg series, 
the so-called forbidden (/) and intercombination (i) lines that together with the 
resonance (r) Is - 2p line form the (ra= 2-1; J = 1 - 0) He-like triplet. For 
0+6, these lines fall at 22.097, 21.804, and 21.602 A, respectively. Among these 
three, r has the highest oscillator-strength by far and, thus, is most affected 
by PE. The intensity ratios of the various 0+^ lines to the r-line as a function 
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Figure 2. Left- Calculated intensity ratios of the Lyman series lines 
and RRC (circles) of 0+^ to the Lyct intensity. The NGC 1068 values 
for these ratios, as obtained from the RGS spectrum, are also plotted 
(asterisks), as well as the values for collisional O^^ gas at 100 cV (hor- 
izontal bars). Right- The NGC 1068 spectrum (dotted line) obtained 
by the RGS spectrometer on board XMM-Newton in the region of the 
0+^ lines and RRC compared with the present model (solid line) at 
10^^ cm-2. 



of column density are plotted in the left panel of Fig. 2. An effect similar to 
that for the H-like species is seen for the high order lines (/?, 7,5, etc.), where 
their intensity with respect to the r-line increases when the latter saturates, and 
reaching a plateau when they saturate as well. The line ratios for 0+^ measured 
from the NGC 1068 spectrum are also plotted in Fig. 2 (asterisks). The same 
as for C"*"^, the measured O"*"^ ratios are found to be higher than the collisional 
ratios, reducing the likelihood of hot collisional gas. The measured ratios of 
Q-l-6 i];npiy a, column density of ~ 3 x 10^^ cm~^. The 0+^ model at this column 
density is plotted versus the data in the right panel of Fig. 2 and, evidently, the 
agreement is fairly good. 



5. Conclusions 



We have shown that a physically consistent model for reprocessing of the nuclear 

continuum and subsequent emission from Seyfert 2 galaxies, which includes all 
of the important atomic processes and particularly PE, reproduces the high- 
resolution RGS X-ray spectrum of NGC 1068 very well. A more complete anal- 
ysis of this spectrum is soon to appear in Kinkhabwala ct al. (2001). Moreover, 
the traditional models, which did not include PE, clearly produce erroneous 
spectra and line ratios. This, in turn, could possibly lead to wrong astrophysical 
interpretation, such as the invoking of hot collisional gas, which to the best of 
our understanding is not evident from the X-rays emitted by Seyfert 2 galaxies. 
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